Engineering artificial graphene in a two-dimensional electron gas 
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At low energy, electrons in doped graphene sheets behave like massless Dirac fermions with a Fermi 
velocity which does not depend on carrier density. Here we show that modulating a two-dimensional 
electron gas with a long-wavelength periodic potential with honeycomb symmetry can lead to the 
creation of isolated massless Dirac points with tunable Fermi velocity We provide detailed theoret- 
ical estimates to realize such artificial graphene-like system and discuss an experimental realization 
in a modulation-doped GaAs quantum well. Ultra high-mobility electrons with linearly-dispersing 
bands might open new venues for the studies of Dirac-fermion physics in semiconductors. 



Graphene is a one- atom-thick two-dimensional (2D) 
electron system composed of Carbon atoms on a hon- 
eycomb lattic^. The lattice has two inequivalent sites 
in the unit cell that are analogous to the two spin 
orientations of a spin-1/2 particle. This observation 
opens the way to an elegant description of electrons in 
graphene as particles endowed with a pseudospin degree- 
of- freedom^. At low energy, electrons in graphene are 
described by a 2D massless Dirac fermion (MDF) Hamil- 
tonian, TYd = vycf • where is the bare Fermi ve- 
locity, which does not depend on carrier density, p is 
the 2D momentum measured from the corners of the 
Brillouin zone, and a is the pseudospin operator con- 
structed with two Pauli matrices {cr*,z = x,?/}, which 
act on the sublattice pseudospin degree-of-freedom. It 
follows that the energy eigenstates are chiral^ i.e. for a 
given p have pseudospins oriented either parallel (con- 
duction band) or antiparallel (valence band) to p. The 
Dirac-like wave equation and the chirality of its eigen- 
states have a number of very intriguing implication^ 
Clearly it would be highly desirable to have other ma- 
terials in which low-energy quasiparticles have Dirac-like 
spectrum and a pseudospin degree-of-freedom. One can- 
didate is represented by HgTe/Hg(Cd)Te quantum wells 
where MDFs are predicted to arise at a critical quantum- 
well thickness^. More recently. Park and Louie^ have 
proposed that MDFs can arise in any 2D electron gas 
(2DEG) if appropriately nanopatterned. 

Here we present an independent approach to the real- 
ization of "artificial graphene" in a nanopatterned 2DEG. 
We provide theoretical evidence for the occurrence of 
linearly-dispersing energy bands in an artificially engi- 
neered honeycomb lattice and we demonstrate a remark- 
able dependence of the Fermi velocity on the strength of 
the external potential in this system. We also define the 
conditions that the external periodic potential and the 



electron density must satisfy in order to achieve artificial 
MDFs. Finally we present the photoluminescence (PL) 
of the 2DEG confined in a high-mobility modulation- 
doped GaAs/AlGaAs quantum well where a nanopat- 
terning with honeycomb symmetry is achieved by dry 
etching. We believe that the development of patterned 
2DEGs with tunable parameters will offer unprecedented 
opportunities to study fundamental interactions of MDFs 
in high-moblity semiconductor structures. 

We start our analysis by considering a 2DEG consist- 
ing of electrons with band mass mb = 0.067 m (m is 
the bare electron mass in vacuum) confined in a thin 
quantum well created in a GaAs/AlGaAs semiconductor 
quantum well. The 2DEG is subjected to a periodic ex- 
ternal potential lext(^) with honeycomb structure. For 
the numerical calculations we have used a 2D muffin-tin 
potential which is zero everywhere but in disks of radius 
r, where it takes the constant value Vq. The center-to- 
center distance between the disks is a, and the lattice 
constant is ao = a/3 a (note that in Ref. |3] they use a 
triangular rather than a honeycomb lattice). 

Because the typical values of ao are much larger than 
the GaAs lattice constant, the external periodic potential 
can be viewed as a long- wavelength superlattice, which 
creates minibands. These are found by solving the secu- 
lar equatiorl^ 
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Here V{G) are the Fourier components of the external 
potential, 



Kxt(r) = ^nG')exp(iG-r) 



(2) 



G = Igi + pg2 are the reciprocal lattice vectors (RLVs) 
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FIG. 1: (Color online) Top left panel: energy minibands (in 
meV and shifted upwards by 0.5 meV) for a repulsive po- 
tential with Vb = +1.0 meV [the points F, M, and K are 
defined in the inset to Fig. ^p)] . Dirac points are formed at 
the K point. The dashed lines are the approximate slopes of 
the bands near the Dirac points. Top right and bottom left: 
same as in the top left panel but for Vb = —0.125 meV and 
Vb = —0.8 meV. Bottom right: the absolute value of the 
threshold potential Vb^^^ (in meV) as a function of r/a, for 
a — 150 nm. 



with I and p integers, and gi^g2 are primitive RLVs, 
gi = 27r(l, V3)/(3a) and g2 = 27r(l, - V3)/(3a). 

In Fig. [l] we plot the calculated energy minibands for 
a muffin-tin potential with a = 150 nm, r = 52.5 nm, 
and for three different values of Vq: Vq = -\-1.0 meV, 
—0.125 meV, and —0.8 meV. As dictated by symmetry 
and group theoTj^, these minibands are characterized by 
the existence of two-fold degenerate points at the corners 
of the Brillouin zone. It is easy to show that states with 
k close to these points are effectively described by a two- 
component MDF Hamiltonian TYd, with a Fermi velocity 
v-p that depends on mb, a, Vb, and on r (see below). 
Comparing the different minibands shown in Fig. [l] we 
clearly see that in the case of a repulsive muffin-tin po- 
tential even when the Fermi level lies at the Dirac point 
(dash-dotted lines in Fig. [T]), other type of states are 
present at the same energy. Ideally, similarly to what 
happens in graphene, one would like to be left only with 
isolated Dirac points at the Fermi level, i.e. with a gap 
in the bulk of the Brillouin zone. Within the muffin-tin 
model we have used, to create such gap we need an at- 
tractive potential, whose strength has to be stronger than 
a certain threshold Vq^^. The absolute value of Vq^^ is 
plotted in Fig. [l] (bottom right panel) as a function of 
the ratio r/a, for a = 150 nm. For a geometry with 
r/a ^ 0.35 the threshold potential is —0.18 meV. 

Note that in the regime where isolated Dirac points 
exist, the band structure of the nanopatterned 2DEG 
at sufficiently low energies consists of manifolds of mini- 



bands separated by minigaps [see e.g. the bottom left 
panel in Fig. [l] . 

Vo, however, cannot be too large in absolute value. In- 
deed, when the local potential Vq is attractive and too 
strong, it can lead to the formation of bound states. In 
this regime any small imperfection in the periodic struc- 
ture, which is experimentally unavoidable, could lead to 
a dramatic change in the character of the states, yield- 
ing complete localization. Transport would occur mainly 
via variable-range hopping, a regime which we want to 
avoid. For a single disk we have estimated this threshold 
potential for the formation of bound states to be 
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meV, 
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where /3 ~ 3.832 is the first zero of the Bessel func- 
tion Ji(x). For r/a ~ 0.35 this localization threshold 
is roughly 3 meV, i.e. twenty times larger (in absolute 
value) than the threshold |Vo^^| necessary to create a gap 
in the bulk of the Brillouin zone. Thus there is a precise 
window of values of |Vo|, which depends on mb and on 
the geometrical parameters a and r, which are suitable 
to create isolated Dirac points in the single-particle band 
structure of the 2DEG. 

When a gap exists in the bulk of the Brillouin zone, 
the condition to have the Fermi level exactly at the Dirac 
point is equivalent to the requirement of having just one 
band filled. The electron density no necessary to fill one 
band corresponds exactly to two electrons per unit cell, 
no = 4/(3V^ a^). For structures with a ~ 150 nm no 
is of the order of a few 10^ cm~^, which is a value that 
can be reached experimentall}'^. For electron densities 
10^ cm~^ < n < 10^^ cm~^, the appropriate a changes 
from 277 nm to 88 nm. 

To explore the applicability of these ideas to real sys- 
tems we realized the periodic external potential on a 
sample containing a 2DEG in a 25 nm wide, one-side 
modulation-doped Alo.iGao.gAs/GaAs quantum well. 
The 2DEG, positioned 170 nm underneath the surface 
(the doping layer is at 110 nm), has measured low- 
temperature electron density Uq = 1.1 x 10^^ cm~^ and 
mobility of 2.7 x 10^ cm^/(Vs). The external modula- 
tion of the 2DEG is achieved following the procedure 
described in Ref. ^ based on e-beam nanolithography 
and inductive coupled plasma reactive ion shallow etch- 
ing (80 nm below the surface). 

Figs. |2|l) and ^) show SEM images of the nanopat- 
terned 2DEG. Other schemes with metallic gates can also 
be explored^. The experimental values of the parameters 
are^ a ~ 150 nm and r ~ 50 nm, similar to those used for 
the calculations. Vq is not known experimentally. How- 
ever, one can engineer systems in which the potential felt 
by the electrons at the corners of the hexagonal cells is 
either repulsive (Vq > 0) or, as in Fig.[2|L) and c), attrac- 
tive (Vb < 0). 

Since Uq > no, electrons occupy not only the bands 
with isolated Dirac points but also minibands at higher 
energy, as shown in Fig. ^), where the density-of-states 
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FIG. 2: (Color online) Panels a) and c) SEM images of 
the nanopatterned modulation-doped GaAs/AlGaAs sam- 
ple. Panel b) Calculated density-of-states (DOS), in units 
of eV~^ X nm~^, as a function of energy for a muffin-tin po- 
tential with Vb — —2.8 meV, a — 150 nm, and r — 52.5 nm. 
The first structure in the DOS is fixed at the energy of the 
main photoluminescence (PL) peak [see panel d)]. The ver- 
tical dashed line denotes the Fermi level associated with the 
nominal density Ue = 1.1 x 10^^ cm~^. The insets show (i) 
the Brillouin zone corresponding to the external periodic po- 
tential represented in panel a) and (ii) a zoom of the DOS 
corresponding to the first structure in the main panel. The 
V-shaped DOS characteristic of MDFs is evident. Panel d) 
Low-temperature (T = 2 K) PL of the sample before (red 
dotted curve) and after (black solid curve) the processing. 



(DOS) and the nominal Fermi level for Vb = —2.8 meV 
are shown. Despite the large doping, the impact of the 
nanopatterning and the formation of minibands clearly 
manifest in the PL spectrum. Fig. [2]i) shows represen- 
tative PL spectra at 2 K both of the unprocessed (red 
dotted curve) and processed (black solid curve) samples. 
In the unprocessed 2DEG case, the PL shape is deter- 
mined by the density-of-states of the free electrons and 
equilibrium occupation factors of the 2DEG and pho- 
toexcited holes^ leading to an estimated electron den- 
sity in agreement with the transport results. The pro- 
cessed sample PL, on the contrary, displays a remarkable 
change with the appearance of sharp structures on the 
high energy side and an overall reduction of its linewidth, 
which are remarkably consistent with the modification 
of the conventional constant-in-energy DOS as shown 
in Fig. ^p)^ provided that Vq is chosen appropriately^^ 
(at Vo = —2.8 meV in this case). The overall smaller 
linewidth of the PL suggests a reduction of the average 
electron density due to the impact of the etching process. 

We would like now to comment on the magnitude 
and tunability of the Fermi velocity of these artificially- 
induced MDFs. Using first-order perturbation theory it 
is possible to calculate analytically the slope of the bands 
close to the K point: in agreement with Ref. 3 we find 




FIG. 3: (Color online) Left panel: color plot of the effective 
Fermi velocity vf [in units of the nearly-free Fermi velocity 
v^^"*] as a function of r/a and |Vb| (in meV), for a = 150 nm. 
Right: vf/v^^"^ as a function of |Vb| for r = 0.35 (circles) and 
T ja — 0.45 (squares), again for a — 150 nm. 

v^"^^ — 27r h / {3 V^mi^a)^ which is {to this order of pertur- 
bation theory) independent of the strength and the sign of 
the potential Vq and is exactly one-half the velocity of a 
free electron of wave- vector K in the absence of the mod- 
ulation. This "nearly free" result applies only for Vb ^ 0. 
In passing, we note that using the value of the bare elec- 
tron mass in vacuum mb = m, and a = 1.42 A (which 
is the Carbon-Carbon distance in graphene) one gets 

v^^"^ ~ 0.98 X 10^ m/s, which is surprisingly close to the 
Fermi velocity of electrons in graphene. Using instead 
the value of the band mass in GaAs, mb ~ 0.067 m, and 
a = 150 nm, we find v^^^ ~ 1.4 x 10^ m/s, roughly an 
order of magnitude smaller than the value of the Fermi 
velocity in graphene. When Vq is beyond the regime of 
applicability of first-order perturbation theory, the gen- 
eral formula for the Fermi velocity [in units of v^^^] is 




where wo{G) are the Fourier components of the peri- 
odic part ofw = 2"^/^(^| - ^|). Here {^a(^), cr =t, j} 
are the two degenerate eigenfunctions of the Hamiltonian 
H = p^/(2mb) + V^xt(^) at the K point, chosen to be a 
basis of the 2D representation of the little group pre- 
sented in Ref. 5 , We have calculated Eq. Q numerically 
and the results are summarized in Fig. |3] We see that a 
finite value of Vb (this plot concentrates only on Vo < 0) , 
away from the regime of applicability of perturbation the- 
ory, tends to reduce the Fermi velocity with respect to 
'^p"^^ . This effect is much stronger here than in Ref. [3j 
we attribute the difference to the different lattice struc- 
ture (honeycomb vs triangular) and to the much smaller 
range of values of Vq considered in Ref. |3l 

Artificially-induced MDFs in 2DEGs confined in high- 
mobility semiconductor heterostructures could offer sev- 
eral advantages over graphene. One is clearly related 
to the very high purity of the former systems. Even 
though the dominant scattering mechanisms in graphene 
are not yet fully understood, it seems that charged impu- 
rities trapped close to (or on) the graphene plane play a 
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very important role in limiting graphene's mobilitjff^ll^, 
partly obscuring intrinsic properties of MDFs. The mass- 
less Dirac fluid at low den sities is indeed a highly inhomo- 
geneous system^^^EIIIIIll^ Although attempts to achieve 
high-mobilities in graphene systems^^ seem to of- 

fer promising prospectives, the possibility of creating ar- 
tiflcially MDFs in high-quality 2DEGs with mobilities 
that can exceed 10^ cm^/ (Vs) even at low densities would 
represent a very exciting alternative rout^^. The real- 
ization of artiflcial graphene thus would pave the way for 
the experimental observation of several predictions made 
for MDFs, such us a universal minimum conductivit}'^ 
o'min = ^e^ /{T^h) and unusual electron-electron interac- 
tion physicp™l|2i. 

Cyclotron-resonance in artiflcial graphene should ex- 
hibit a ^/B dependence and quite large gaps. Indeed, 
the MDF cyclotron frequency i^ ujf^^ V^^f/^b, 
where = \/hc/{eB) ^ 257 A/^^^/Tesla is the mag- 
netic length, while in a standard 2DEG it is instead 
^2D EG ^ h/{m\^l\). In graphene vy ~ 10^ m/s and 
thus, at a fleld 5 = 10 T, uf^^ - 1328 K. In GaAs, at 
the same fleld, uo'^^ ~ 63 K. In artiflcial graphene cre- 
ated in GaAs, as we have seen, the Fermi velocity is tun- 
able: if we use its nearly- free value, ~ 1-4 x 10^ m/s, 
we get uof^^ - 186 K at 10 T, which is roughly one order 
of magnitude smaller than in graphene but still one order 
of magnitude larger than in a standard 2DEG. More im- 
portantly, since the fractional quantum Hall effect is eas- 
ily observed in 2DEGs, artiflcial graphene could also be 
a very useful playground to understand electron-electron 
interaction effects in MDF systems in the presence of high 
magnetic flelds. The large area of the superlattice unit 
cell offers also the possibility to achieve commensurabil- 



it}/^^ between an external magnetic flux and the quantum 
of flux at quite small values of the external magnetic fleld. 

The realization of MDFs in conventional semiconduc- 
tors opens the interesting scenario related to the im- 
pact of spin-orbit coupling particularly if one uses InAs- 
based materials where spin-orbit coupling and effective g- 
factors are large. Finally we would like to mention that 
if the nanopatterning technique discussed above is car- 
ried out on areas with extension much smaller than that 
illustrated in Fig. |2j one can in principle create artiflcial- 
graphene ribbons with perfect zigzag or armchair edges in 
which conflnement and boundary conditions could play 
a very important role. 

In summary, we have shown that under suitable condi- 
tions systems of massless Dirac fermions can be created 
in conventional 2DEGs conflned in semiconductor quan- 
tum wells by creating a superlattice with honeycomb ge- 
ometry by nanopatterning. The existence of artiflcially- 
induced masslessness could be demonstrated through the 
observation of the half-integer quantum Hall effect^ or by 
studying the PL or neutral collective excitations in the 
low-density regime. Such artiflcial graphene structures 
embedded in semiconductors could open novel routes for 
studies of electron interactions in low-dimensional sys- 
tems. 
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